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Abstract. The electronic structure of Bernal-stacked graphite subject to a tilted 
magnetic field is studied theoretically. The minimal nearest-neighbor tight-binding 
model with the Peierls substitution is employed to describe the structure of Landau 
levels. We show that while the orbital effect of the in-plane component of the magnetic 
field is negligible for massive Dirac fermions in the vicinity of a the K point of the 
graphite Brillouin zone, at the H point it leads to the experimentally observable 
splitting of Landau levels, which grows approximately linearly with the in-plane field 
intensity. 
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1. Introduction 

The recent attention paid to graphene monolayers has been motivated by their unusual 
two-dimensional (2D) Dirac energy spectrum of electrons. In Bernal-stacked graphene 
multilayer composed of weakly coupled graphene sheets, the interlayer interaction 
converts the 2D electron energy spectrum of graphene into the three-dimensional (3D) 
spectrum of graphite. The electronic structure of 3D graphite subject to magnetic fields 
perpendicular to x — y planes of graphene layers were extensively studied a long time 
ago, see, e.g., Refs. [B [21 S H El [6] . 

The application of the tilted magnetic field B = {0,By,Bz) is a standard method 
used to distinguish between 2D and 3D electron systems, as in 3D systems the orbital 
effect of the in-plane magnetic field component should be observable. We will study 
this problem theoretically using a simple tight-binding quantum mechanical model of 
the graphite electron structure. 

Various approaches were employed previously to study the influence of the tilted 
magnetic fields. 

The Fermi surfaces of metal single crystals were investigated by measurements of the 
de Haas-van Alphen effect in tilted magnetic fields. The interpretation of experiments 
relies on the quasiclassical Onsager-Lifshitz quantization rule [3 [8], the Fermi surface 
is reconstructed from the periods of magneto-oscillations which are proportional to 
angular-dependent extremal-cross-sections perpendicular to the direction of the tilted 
magnetic field. 

In semiconductor superlattices the quasiclassical interpretation of data measured 
in tilted magnetic fields fails, as reported in Refs. [9l [TOl [HI [T21 [131 [H]. In these papers 
the observed quantum effects are attributed to the shift of centers of fc-space orbits 
in neighboring quantum wells by \e\By/d in the /c^i-direction, where d is the distance 
between quantum wells. In real space this means that the in-plane magnetic-field length 
iy = ^Jh/\e\By should become comparable with d to reach the visible effect [15] . 

Besides semiconductor superlattices, other layered materials with much shorter 
interlayer distances were also investigated in tilted magnetic fields. Different versions of 
angular magnetoresistance oscillations (AMRO) were studied both experimentally and 
theoretically in low-dimensional quasi-2D and quasi-one-dimensional organic conductors 
(see, e.g., Refs. [T6l[I71[T8l[l9l[2Ql[2Il[221[23l[21] and references therein), and also in 
intercalated graphite [251 [26]. On the theory side, the high Landau level (LL) filing 
factors and weak interlayer interaction were considered in iterpretation the data. 

In pristine graphite this problem has been touched on by two recent theoretical 
articles. 

The graphene multilayer energy spectrum in magnetic fields parallel to the layers 
was described quantum-mechanically in Ref. [27], as the standard theory of AMRO in 
tilted magnetic fields was not applicable due to the relatively strong interlayer interaction 
(in comparison with the intercalated graphite) between graphene sheets. 

The LLs in the bilayer graphene in magnetic fields of arbitrary orientations were 
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calculated analytically in Ref. [25] . 

Both papers conclude that a very strong in-plane field component is necessary to 
induce an observable effect on the electronic structure. Indeed, to reach iy comparable 
with the distance between graphene layers in graphite, the magnetic field By = 5865 T 
would be necessary. 

In this paper we make use of the specific features of the LL structure in two 
nonequivalent neighboring graphene sheets in graphite, and show that at the H point of 
the graphite hexagonal Brillouin zone the application of the tilted magnetic field leads 
to experimentally observable splitting of LLs of the order of several meV. 

2. Model 

Bulk graphite is composed of periodically repeated graphene bilayers formed by two 
nonequivalent Bernal-stacked graphene sheets, as shown in Fig. [H There are two 
sublattices, A and B, on each sheet and, therefore, four atoms in a unit cell. The 
distance between the nearest atoms A and B in a single layer is 1.42 A, the interlayer 
distance between nearest atoms A is d = 3.35 A. 




Figure 1. (Color online) The lattice structure of graphite. The unit cell is a green 
parallelepiped. 

To describe the graphite band structure, we employ the minimal nearest-neighbor 
tight-binding model, introduced by Koshino et al. in Ref. [29]. This model is reduced 
Slonczewski, Weiss and McClure (SWM) model. The tight-binding Hamiltonian "H 
includes only two instead of seven tight-binding parameters, the intralayer interaction 
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7o = 3.16 eV between the nearest atoms A and B in the plane, and the interlayer 
interaction t = 0.39 eV between the nearest atoms A out of plane. 

While the reduced SWM model is not appropriate, e.g., for a Fermi surface 
description, this model has been successfully applied in the theoretical papers [271 128] . 
and used to describe recent magneto-optical measurements on graphite in Refs. [3ll 
[32l [30l [33| [M] . It has been shown in Refs. [301 [331 IM] that the transitions between 
Landau levels originating from the K and H points of the graphite Brillouine zone can 
be understood within a simple picture of an effective bilayer with a coupling strength 
enhanced twice in comparison to a true graphene bilayer, 2t, (which is definitely the 
effect of a superlattice) and an effective graphene monolayer. These solid arguments are 
in agreement with the theoretical model based on the reduced SWM model we develop 
in our manuscript. 

In this model the wave functions are expressed via four orthogonal components ipf, 
ipf , "^i+i' '^f+i^ which are, in zero magnetic field, Bloch sums of atomic wave functions 
over the lattice sites of sub lattices A and B in individual layers j. 

The continuum approximation is used in the vicinity of the H — K — H axis of the 
graphite hexagonal Brillouin zone, for small k = {k^, ky) measured from the axis. Then 
the electron wave length is larger than the distance between atoms, and the non-zero 
matrix elements of "H can be written as 

H^^ = hvF{k^ + iky), (1) 
U''^ = hvF{k:,-iky), (2) 
H^^ = t. (3) 

The Fermi velocity, vp, is defined by Uvp = \/3a'jo/2, and will be used as an intralayer 
parameter instead of 70 in the subsequent consideration. 

The effect of the arbitrary oriented magnetic field, B = (0, By, B^), can be 
conveniently introduced into the zero-field Hamiltonian by the Peierls substitution. If we 
choose the vector potential in the Landau form A = {ByZ — B^y, 0, 0), the substitution 
will read 

hkj; hk^ - \e\B^y + \e\Byjd, (4) 

where an integer number j indicates the graphite layer number. Consequently, the 
matrix elements, Ti^^ and H^^, become layer dependent in a tilted magnetic field, 

Hf^ = vpikx - \e\Bzy + \e\Byjd + iky) = vpHj, (5) 
nf^ = VF{k, - \e\B,y + \e\Byjd - iky) = vpU*. (6) 

Making use of the above approximations, the Schrodinger equation involving all 
layers j, leads to the following system of equations 

VFU*^f - E^f + t^f^,+t ^f^, = 0, (7) 



-E^pf + VFTl.*i,f = 0, (8) 
VFll*+,^f^r - E^f^, +t^f + t = 0, (9) 
-E'i(:f^,+VFll,+,ijf^, = 0. (10) 
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Note that the structure of H allows to express the function ip^ via the function ipj^ from 



the same layer, and we are thus left with two interlayer equations for ip^ and ip^ 



(4n,n* 



E')^pf + tE{^f_, + i;f^,) = 



(4n*+in,+i - E')^pf^, + tE{^f + ^f^^) 
3. Zero-field case 



0. 



(11) 
(12) 



It follows from the condition of periodicity in the ^-direction that ipfj^i and ipf can be 
written as 



where (p^ and (p2 denote the 2D wave functions in two non-equivalent layers of the 
graphite unite cell, and is restricted to the first Brillouin zone, —tt/2 < k^d <ti/2. 
For i? = Eqs. (fTTl [T2l) are transformed to 

[h^vlk'^ - E^)(P^ + 2tE cos{k,d)(pf = 0, (14) 



JkzdjiA 



(13) 



{h'v'pe - E' 



'4 



2t E cos{k,d)(p2 = 0, 



and from here the four eigenvalues are obtained 



E^'^ = ±r± + h^ik^, 

where 

T = t cos{kzd) 

denotes the /c^-dependent coupling of two graphene sheets. 



(15) 
(16) 
(17) 



4. Perpendicular magnetic field 

In the perpendicular magnetic field. By = 0, the system remains periodic in the z- 
direction, and Eq. f lT5]) is still valid. To find a proper form of T-L^^ and "H^"^ we introduce 
the perpendicular magnetic field length, = h/{\e\Bz), the centre of a cyclotron orbit. 
Ho = i^kx, the dimensionless variable, rj = {y — yo)/iz, and the perpendicular-magnetic- 
field-dependent parameter, B = 2h\e\BzVp. Then, in our notation. 



n 



AB 



H 



BA 



B 

2 

B_ 

2 



d_ ^ 
dr] 



'd_ 
dr] 



and a being raising and lowering operators, respectively. 
With help of these expressions, Eqs. (fTTj) and f|T2l) can be written as 



'B 




2 




f( 









E' 



E' 



+ 2r^0^ = 0, 



'^^ + 2rE<p^ = 0. 



(18) 
(19) 



(20) 
(21) 
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It is obvious from these equations that (pf and 0^ are closely related to the eigenf unctions 
of the harmonic oscillator, ipniv)- Assuming 



1 



n'=0 



n'=0 



and having in mind that 
/ 92 



we get 



B{n + 1) - E^] A2,n + '2EY^ T^x^ix 

n'=0 



Bn - 

where T„^„/ is defined by 



n'=0 



0, 



(22) 

(23) 

(24) 
(25) 



/ + CXD 
Vn{'n)'Pn'{ri)dr] = T6n,n'. (26) 
~oo 

It follows from fl26|l that only the LLs with the same quantum numbers n (but with the 
different energies) are coupled and we arrive to 

'B{n + 1) - E^] A2,n + 2Tn,nEAi, 



0, 



Bn - E^] Ai„ + 2Tn,nEA2, 



0. 



Solving Eqs. ( ITTt [25]) yields the eigenenergies 



± 



(27) 
(28) 

(29) 



which are presented in Fig. O 

As the densities of states of the above Landau subbands have singularities at the 
points K and H, kzd = and k^d = 7i/2, respectively, we concentrate on the field 
dependence of levels corresponding to these points. 

At the K point the eigenenergies are given by 



E. 



± 



2t^ + B 



1\ / 1 B"^ 

n + -j±J4t^ + 4t^B{n + -) + -., 



(30) 



i.e., the subbands are equivalent to those of a graphene bilayer with the coupling constant 
doubled: 2t instead of t. 

At the H point the eigenenergies read 



Et^ = ±\lB[n + -jT ^ 



B 



(31) 



Due to the effectively vanishing inter-layer coupling the spectrum corresponds to the 
Dirac fermions. The coefficients Ai^n and A2,n in Eqs. ( 1271 1281) are equal to 1 and this 
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0.7854 1.5708 2.3562 3.1416 

Figure 2. (Color online) Landau subbands of graphite, E^'^ and E^-^ (denoted as 
and £'^„ at the H point) subject to perpendicular magnetic field i?2 = 5 T, as a 
function of kz along the K-H-K pass of the Brillouin zone. The fc^ dependence of LLs 
in the vicinity of the H point restricted by the rectangular box will be compared with 
those subject to the tilted magnetic field, as shown in Fig. |6l 

implies that the corresponding wave functions are locahzed either in the layer 1 or in 
the layer 2. To emphasize it, we will write the eigeneneries as 

= = ± VB^, (32) 

E±„ = = ±^Bin + l). (33) 

The energy spectrum of LLs at the H point is presented in Fig. [31 

We also denote the wave functions 4>^2) and 0f(2) as |?t.)^2) and |?^)f(2) to stress that 

they are the envelope wave functions of atomic wave functions A and B in the layers 1 

and 2, as shown in TABLE [H 

Let us mention that E^^ = E'f^^^, i.e., we are left with pairs of degenerated LLs 

with different quantum numbers n but the same eigenenergies, with the wave functions 

localized in two different layers. 

5. Tilted magnetic field 

While the previous two paragraphs summarized the already published theories devoted 
to i? = and B = (0, 0, B^), here we present new results for B = (0, By, B^). 
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B,(T) 

Figure 3. (Color online) LLs of graphite at the H point of the Brillouin zone, -Ej^, 
and in the perpendicular magnetic field, Bz- 



Table 1. Energies and wave functions at the graphite H point in perpendicular 
magnetic fields 



Energy 


{<Pt < 


/>f, 0f} 


Et = VB^ 


{\n)t- 


|n-l)f,0 ,0} 


= -^/Bn 


{\n)t 


|n-l)f,0 ,0} 


Et = ^B{n + l) 


{0, o,h 


n)l-\n + l)^] 


E:^ = -^B{n + l) 


{0, 0,1 


n)t |n + l)f} 
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5 10 15 



Figure 4. (Color online) The dimensionless parameter rjd as a function of the 
perpendicular component of the magnetic field, B^, calculated for various tilt angles 
of the magnetic field B. 



In tilted magnetic fields the off-diagonal matrix elements, "Hf^ and "Hf^, given by 
Eqs. (El E]) remain layer dependent and take the form 



H 



AB 



BA 




(34) 



(35) 



(36) 



where the small dimensionless parameter 
By d 

means the shift of the cyclotron orbit center in the j-layer due to the in-plane component 
of the magnetic field, By. Fig. H] illustrates the B^ dependence of r]^ for various tilt angles 
of the magnetic field. 

Note that graphite subject to tilted magnetic fields is no longer periodic in the z- 
direction, but becomes periodic in the direction of the tilted magnetic field. To take into 
account the shift of the cyclotron orbits, we apply the approach developed in Ref. [35] 
for semiconductor superlattices. Accordingly, the Eqs. (ITTl) and ( IT2l) are modified to 





Q2 










5( 





+ iv- U - 1)%)^ 



0, 



^^1 + 



+tE{i;f + i;f_,,) 



0. 



(37) 



(3J 



Tight-binding description of Landau levels of graphite in tilted magnetic fields 10 



The new periodicity implies that ip^^ and "ipf^i can be written as 



= e^^-''^^ iv + m) , = [r/ + (j + l)r/,] , (39) 

where again —tx/2 < k^d < 7r/2. Here the wave functions 0^ (r/ + jrjd) and 0^ (77 + jrjd) 
are associated with cyclotron orbits in two layers. Introducing the shift operator by 



0(^ + %) = e*'"'^0(^), K = -id/dT], 



(40) 



and employing the ^-representation, two interlayer Eqs. ( 1371) and (l38l) can be given the 
form similar to Eqs. (1201 [2T]) for the perpendicular magnetic field 



2" 
B 
2 



1 -E^ 



+ 2r^0f = 0, 



b^ + 2r^0^ = 0, 



(41) 
(42) 



but with the new coupling T which depends, due to new periodicity, not only on kz but 
also on both components and By of the arbitrary oriented magnetic field via r/^, 

T(k) = t cos {KTjd + kzd). (43) 

The Eqs. (14T!) and (142|) represent the main result of this paper and in the following 
we will discuss the possible methods of their solutions. 

Expressing again 0^ and with a help of (EH) and ([25]) we arrive to 



B{n + 1) - 



n.'=0 



0, 



Bn-E' 



^1, 



2-E ^ Tn,n'A2,n' 
n'=0 



0, 



where 



(Pn{K)T{K)iPn'{l^)dn. 



(44) 
(45) 

(46) 



The integrals (1461) can be evaluated analytically and expressed via the generalised 



Laguerre polynomials (see, e.g., Ref. [36] )• The Eqs. (jH]) and (H5]) define the matrix 
which should be diagonalized. The nonzero coupling of LLs with different n allows us 
to conclude that the degeneracy of LLs at the H point will be removed, and the LLs 
with different n avoid to cross. Then the standard approach is to solve the secular 
equation numerically, the minor complication being that the matrix elements depend 
on the energy. 

Here we prefer to obtain analytic results by the lowest order perturbation theory 
application and treating 77^ as a small parameter. 

A simple trigonometric relation and a series expansion restricted to terms linear in 
r]d imply 

T^k) = t cos{kzd) cos {K,r]d) — t sin (kzd) sin (^77^) ^ 

^ t cos{kzd) —t sin {kzd)ridH- ■■ , (47) 
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and 



(Pn{K)K(pn'{K)dK. 



(4^ 



Tn,n' = t COs{k^d)6n,n' " % sill (k^d) 

The integrals in Eq. ( 148|) can be easily evaluated using the relation (see, e.g., Ref. [37] ) 

(49) 



ln + 1 



V 2 "" ^ ' V 2 

Let us pay attention to the field dependence of LLs at the most interesting K and 
H points of the graphite Brillouin zone. 

At the K point fc^rf = 0, and, consequently, cos{kzd) = 1, sm{kzd) = 0. Then Tny 
reduces to 



Tn. 



tSn 



(50) 



This coupling corresponds to the effective bilayer subject to the perpendicular field 
discussed above. Thus, we have found that corrections induced by By are very small 
and of the order 77^. This is in agreement with conclusions presented in Ref. |28] . 

The field dependence at the H point, k^d = ±vr/2, is more interesting. In 
perpendicular magnetic fields the coupling between layers disappears, and, as presented 
above, we obtained the LLs corresponding to graphene Dirac fermions, namely = 
±\/Bn for the first layer, and E^^ 



±^B{n + 1) for the second layer, n = 0, 1, 2, ■ ■ -. 
In the magnetic field of an arbitrary direction the interlayer interaction is not 
reduced to zero, but remains finite. The non-zero matrix elements Tn,n' can be written 

as 



'Tn,n-\ 



l,n 



trjd 



In 



Tn,n—1 "Tn- 



l,n 



In 



(51) 



The small perturbation rj^ couples the states with \n + 1)^ and |n — 1)^. 
Among them the states |n)^ and |n + 1)^ are degenerated, i.e., they belong to the 



same unpertubed eigenvalues ±^^jB{n + l). Consequently, at least the lowest order 
perturbation approach suitable to remove the degeneracy must be applied, which yields 
equations 

'B{n + 1) - E^] A2,n + 2Efn,n+lAi,n+l = 0, (52) 
'B{n + 1) - E^] + 2Efn+l,nA2,n = 0. (53) 

The secular equation derived from Eqs. (!52l [53]) reads 



B{n + 1)-E'\ -AE'T^^^, 
and from here we get the four eigenenergies 

± 



0, 



(54) 



„+i - ±r„,„+i ± ^JB{n + 1) + r„%+i, n = 0, 1, 2, ■ ■ ■ (55) 

The eigenenergies, E"^ and E^^ steming from Ef^^ remain the same as in the 
perpendicular magnetic field. In that case the degeneracy is not removed. 
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Figure 5. (Color online) Splitted by tilting the magnetic field with ip = 20° LLs of 
graphite at the H point, E^'^, as a function of the perpendicular component of the 
magnetic field, B^. The dotted lines correspond to results of numerical calculations. 



Lifting of LL degeneracy by the tilted magnetic field in LLs with n > is shown in 
Fig. [5l The LL splitting is of the order of several meV, and it grows with the tilt angle, 
i.e., with the in-plane magnetic field component. By. 

The corresponding eigenfunctions calculated with the same level of accuracy are 
presented in TABLE [2j They are mixed from wave functions of both layers with an 
equal weight. 

To test the accuracy of the above approximations we have calculated the eigenvalues 
numerically with the larger basis \n)i, \n + 1)^, \n + 2)^, jri)^, \n + 1)^, \n + 2)^ instead 
of the minimal one \n + At the H point we have found only negligible 

quantitative corrections to the results obtained analytically, as presented in Fig. |5l 
Similarly, including the higher order expansion in rjd does not influence the results for 
the chosen range of angles and magnetic fields. 

The above basis allows to calculate also the dependence in the vicinity of the 
H point defined roughly by the rectangle in Fig. [2j The most interesting feature is the 
development of additional local extrema near the H point, which are more pronounced 
for LLs with higher n. The same is true for minigaps open at the crossing points of LLs, 
as mentioned in the previous paragraphs. The results are shown in Fig. |6l 

In general, our approach must fail for magnetic fields close to the in-plane 
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Table 2. Energies and wave functions in tilted magnetic fields at the graphite H point 



Energy 








{\n)t-\n-l)f,- 


n + |n + 2)f} 


E-- 


\n - l)f , -\n 


+ l)^,-|n + 2)f} 


Et^ 


{\n + l)t-\n)f, 


n)^,-|n + l)f} 


Et 


{\n + l)t\nr,, 


r2)^,|n + l)f} 



0.22 



0.2 s 



0.18 



0.16 
Uj 0.14 



0.12 



0.1 



0.08 - 



1 .5237 



e/' rep = 20°) 



H 



_L 



_L 



1 .5394 



1.5551 



1.5708 

k d 

z 



1 .5865 



1.6022 



n = 5- 
n = 6. 



n = 4' 
n = 5 






1.6179 



Figure 6. (Color online) Graphite Landau subbands, depicted by the rectangular box 
in Fig. [5J as a function of nearby the H point. Dashed curves are graphite LLs, 
and E2n7 in the perpendicular magnetic field = 5 T, = 0°. Solid curves 
are LLs of graphite, En^^ and En^^ , subject to the tilted magnetic field with the 
perpendicular component -B^ = 5 T and the tilt angle = 20°. 



Tight-binding description of Landau levels of graphite in tilted magnetic fields 14 



orientation, By ^ B^, a,s rj^ oo for B^ ^ 0, and the expansion in powers of rj^ 
is no longer acceptable. 

Also the perturbation theory is less appropriate for states with large n, as the 
energy difference between neighboring LLs is smaller then for states with small n and, 
moreover, the interlayer coupling matrix elements, Tn,n±i, increase with ^/n. The limits 
of the numerical approach are not clear at present, but we should have in mind that 
from the experimental point of view the angles with almost in-plane orientation are not 
so interesting due to the mosaic structure of most graphite crystals. 

6. Conclusion 

Based on the simple nearest-neighbor tight-binding quantum mechanical model, we 
presented the calculation method of the band structure of Bernal-stacked graphite 
subject to tilted magnetic fields. We applied the lowest order perturbation theory to 
obtain analytic solutions of the formulated equations, the accuracy of which was later 
checked by the simplified numerical calculation. The special attention has been paid 
to the field dependence of the LLs at the K and H points of the graphite Brillouin 
zone where the density of states exhibits van Hove singularities in the perpendicular 
magnetic field. We have found that at the K point, where the electron structure in the 
perpendicular magnetic field reminds strongly that of the bilayer graphene, the influence 
of the in-plane component of the magnetic field is negligible. On the other hand, at the 
H point, where the electron structure mimics the behavior of the Dirac fermions, the 
application of the tilted magnetic field leads to the splitting of LLs. This splitting is 
of the order of several meV, which is an experimentally observable value, and it grows 
with increasing of the in-plane component of the magnetic field. 
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